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Aerodynamic  drag  and  lift  effects  on  the  inclination  of  the  orbit  of  the  S3-1 
satellite  have  been  used  to  infer  atmospheric  zonal  wind  speeds  at  an  altitude 
of  175  km  and  to  estimate  the  satellite  surface  thermal  accommodation  coef- 
ficient. The  unusually  high  quality  of  the  orbital  and  attitude  data  permitted 
the  7 -month  data  span  to  be  analyzed  in  6 non-overlapping  subintervals, 
thereby  resulting  in  finer  local-time  and  spatial  resolution  than  is  normally 
obtained  in  such  studies.  The  aerodynamic  lift  force  was  found  to  be  quite  — 
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sensitive  to  the  assumed  value  of  the  thermal  accommodation  coefficient, 
a quantity  that  is  neither  well  known  nor  understood.  The  inclination  data 
were  best  fit  (in  a least  squares  sense)  when  the  accommodation  coefficient 
was  greater  than  0.95.  Best  agreement  with  theoretical  models  in  terms 
of  wind  speed  and  direction  was  obtained  when  the  accommodation  coeffi- 
cient was  taken  to  be  nearly  1.0.  In  this  case,  the  inferred  wind  was  ^-400  m/ 
sec  eastward  at  local  times  19-21  hr,  and ^-200  m/sec  westward  at  local 
times  7-9  hr.  The  study  indicates  that  if  accurate  zonal  wind  velocities  are 
to  be  derived  from  satellite  orbital  inclination  changes,  care  must  be  taken 
to  ensure  that  the  shape  and  orientation  of  the  satellite  are  such  that  the  lift 
forces  generated  are  sufficiently  small  that  the  resulting  uncertainty  due  to 
inadequate  knowledge  of  the  accommodation  coefficient  does  not  cause 
significant  error. 
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INTRODUCTION 


The  analysis  of  long-term  changes  in  the  inclination  of  earth- 
orbiting satellites  is  a valuable  tool  for  investigating  the  characteristics 
of  upper-atmospheric  winds.  The  majority  of  this  work  has  been 
carried  out  by  King-Hele  and  his  colleagues  (see  King-Hele,  1972a 
for  a review)  and  has  revealed  an  average  eastward  rotation  of  the 
atmosphere  in  the  altitude  range  200-350  km.  Little  information 
regarding  the  north-south  component  of  the  winds  has  been  gathered 
in  this  way  because  of  the  small  effect  of  this  wind  component  on  long- 
term satellite  inclination  changes;  however,  recent  studies  (King-Hele, 
1976)  have  been  aimed  at  this  aspect  of  the  problem.  Furthermore, 
the  diurnal  variations  of  the  wind  field  have  been  treated  in  this  way 
only  very  recently  (Blum  and  Schuchardt,  1976),  but  this  study  applied 
only  to  exospheric  heights.  The  primary  obstacle  in  such  a study  arises 
from  the  fact  that  the  local  time  at  perigee,  where  the  forces  normal 
to  the  orbit  plane  due  to  winds  are  most  effective,  changes  considerably 
on  the  time  scale  necessary  for  significant  changes  in  inclination  to 
occur;  thus,  the  diurnal  variations  in  winds  are  difficult  to  isolate. 
However,  in  principle,  this  technique  is  complementary  to  incoherent 
scatter  radar  (Evans,  1972;  Amayenc,  1974)  and  Fabry-Perot  inter- 
ferometric (Hernandez  and  Roble,  1976)  determinations  of  thermospheric 


I 
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None  of  the  previous  studies  of  satellite  inclination  changes  has 
taken  explicit  account  of  possible  effects  of  aerodynamic  lift.  This 
is  primarily  because  no  accurate  attitude  information  has  been 
generally  available.  Furthermore,  the  aerodynamic  lift  effects 
on  spherical  or  near-spherical  satellites  should  be  negligible. 

However,  Karr  and  Yen  (1970)  and  Karr  et  al.  (1975)  have  shown  that 
the  effects  of  aerodynamic  lift  on  various  orbital  elements  can  be  consid- 
erable, especially  for  a satellite  which  has  large  planar  surfaces  or  a 
fixed  attitude  relative  to  its  velocity  vector. 

The  purpose  of  the  present  paper  is  to  evaluate  the  effects  of 
upper-atmospheric  winds  and  aerodynamic  lift  on  the  inclination  of  the 
S3-1  satellite.  As  will  become  clear,  the  facts  that  the  orbit  was  nearly 
sun-synchronous  and  that  the  orientation  of  the  satellite  was  almost 
fixed  and  was  closely  monitored  allow  a detailed  study  of  diurnal  variations 
of  upper-atmospheric  winds  and  the  effects  of  aerodynamic  lift  on  the 
inferred  values  of  wind  speeds,  and  an  estimate  of  the  thermal  accommodation 
coefficient. 
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THE  SATELLITE  AND  ITS  ORBIT 


S3-1,  a research  satellite  developed  under  the  U.  S.  Air  Force 
Space  Test  Program,  was  placed  into  near  sun- synchronous  orbit  in 
late  1974.  The  payload  consisted  of  a number  of  instruments  designed 
to  measure  thermospheric  density  and  composition.  On  November  4, 
1974,  characteristics  of  the  orbit  were  as  follows: 


Apogee  height 
Perigee  height 
Eccentricity,  e 
Inclination,  i 

Right  ascension  of  ascending  node,  ft 
Argument  of  perigee,  u: 

Period,  T 


3767  km 
159  km 
0.  217 
96.98  deg 
29.04  deg 
125.43  deg 
126.1  min. 


Figure  1 shows  the  evolution  of  i,  uu  and  T during  the  period 
November  4,  1974  to  May  26,  1975.  Since  the  orbit  is  nearly  polar, 
the  latitude  of  perigee  X can  be  determined  by  inspection  of  u>  in  Fig.  1 
(sin  X = sin  uu  sin  i).  The  local  time  of  perigee  occurred  approximately  in 
the  meridian  plane  containing  the  local  times  8 and  20  hr.  (Further  infor- 
mation regarding  the  local  time  of  perigee  is  given  in  Table  2 below.  ) 

The  orbit  of  the  S3-1  satellite  was  especially  suitable  for  drag 
analysis  because  of  the  unusually  high  quality  and  completeness  of  orbit 
and  attitude  parameters.  Tracking  was  performed  by  the  USAF  S-band 
Space -Ground  Link  System  (SGLS),  which  utilizes  an  up-link  active  tracking 
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signal  having  a pseudo-random  noise  (PRN)  ranging  code  which  is 
echoed  back  to  the  remote  tracking  station  by  a phase-locked  transponder. 
Round-trip  timing  of  the  signals  along  with  Doppler  shift  measurements 
allowed  range  and  range-rate  to  be  derived.  These  quantities  along  with 
the  observed  satellite  azimuth  and  elevation  were  used  to  compute  orbital 
elements.  On-board  earth  horizon  sensor,  sun  angle  sensor,  and  three- 
axis  magnetometers  allowed  determination  of  the  satellite  attitude. 

Magnetic  torque  coils  enabled  the  spin  axis  to  be  maintained  nearly 
normal  to  the  orbit  plane.  The  satellite  was  contacted  about  six  times 
per  day  in  order  to  ensure  adequate  tracking  information.  Orbital  elements 
and  spin  axis  orientation  were  computed  at  least  once  per  day  throughout 
the  satellite  lifetime.  The  inclination  data  has  a root-mean-square  (rms) 

_3 

error  on  the  order  of  10"  deg.  Two  hundred  thirteen  sets  of  orbital 
elements  and  spin  axis  orientation  parameters  (right  ascension  and 
declination)  were  used  in  the  present  study. 

In  addition  to  the  orbital  elements  and  attitude  parameters,  the 
satellite  shape  is  of  importance  in  considering  aerodynamic  effects. 
Photographs  of  the  actual  satellite  and  access  to  an  accurate  full  sized 
mock-up  ensured  adequate  knowledge  of  shape  and  dimensions  for  this 
study. 
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CALCULATION  OF  AERODYNAMIC  FORCES 

The  aerodynamic  forces  exerted  on  a body  traveling  through  a 
gas  in  the  free  molecular  flow  regime  depend  on  the  details  of  the 
interaction  of  the  impinging  gas  particles  with  the  surface  of  the  body. 

A two-parameter  model  (Karr  and  Yen,  1970;  Karr  et  al.  , 1975)  has 
been  used  to  describe  this  interaction.  The  forces  acting  on  a flat  plate 
of  area  S are  illustrated  in  Fig.  2.  As  viewed  from  the  plate,  gas  particles 
with  average  velocity  v^  approach  from  a direction  making  an  angle  x with 
the  surface  normal.  The  velocity  and  angular  distribution  of  the  reflected 
particles  are  characterized  by  two  parameters,  a and  P,  respectively. 

The  velocity  of  the  reflected  particles  is 


v = v 
r ° 


v/1  ' a > 


(1) 


so  that  a resembles  the  usual  thermal  accommodation  coefficient, 
which  appears  in  the  relationship 


2 2 
v - v 

r * 


(1  - 


°T' 


+ vs  aT 


(2) 


The  velocity  v , which  characterizes  the  surface  temperature,  Tg , of 
7 s 

the  plate,  is  given  in  terms  of  k,  Boltzmann's  constant,  and  m,  the 
molecular  weight  of  the  gas  molecules,  by  v 2 = 3kT  / 2m.  The  average 
angle  of  reflection  is 
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r 


x r = x (1  - P). 


and  the  angular  distribution  of  the  reflected  particles  ranges  from  pure 
specular  at  P = 0 to  diffuse  at  P - 1. 

The  aerodynamic  force  on  the  plate  may  be  resolved  into  a drag 
component,  F~,  parallel  to  v , and  a lift  component,  F , perpendicular 
to  v and  in  the  plane  of  v and  the  normal  to  the  surface,  given  by 


Z cos  x - Z J\  - a cos  x cos  [ v - (Z  - P)  x 3 


(4) 


Fl  = I P vf  S 


Z sj  1 - a cos  x sin  [tt  - (Z  - P)  X ] 


where  p is  the  mass  density  of  the  gas.  The  quantities  in  the  braces  of 
Equations  (4)  and  (5)  are  the  so-called  drag  (Cp)  and  lift  (C.jJ  coef- 
ficients, respectively. 

The  S3-1  satellite,  schematically  illustrated  in  Fig.  3 is  approximated 
as  an  assembly  of  eight  flat  plates.  In  orbit  it  was  maintained  with  its 
spin  axis  nominally  perpendicular  to  the  orbit  plane  with  a spin  rate  of 
about  5 rpm.  Equations  (4)  and  (5)  may  be  used  to  calculate  the  lift  and 
drag  forces  averaged  over  one  spin  rotation  of  the  satellite. 

Let  x,  z be  mutually  orthogonal  unit  vectors  defining  a rectangular 
coordinate  system.  The  spin  axis  is  chosen  to  lie  along  the  z axis.  Consider 
a pair  of  identical  flat  plates,  each  of  area  S,  symmetrically  placed  on 
opposite  sides  of  the  z axis.  Let  the  outward  unit  normal  vectors  for 
the  plates,  ri  and  n be  given  by 

”“l  ~3 
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nj  cos  ^ sin  T]  x,  + sin  5 sin  T|  y + cos  7]  jz  , 


<N  _ _ /N  . /N. 

n^  - cos  r sin  t|  x - sin  I?  sin  T)  y + cos  T)  z. 


where  £ and  t|  are  the  azimuthal  and  polar  angles  of  n respectively. 


For  5 = 0 the  geometry  is  as  illustrated  in  Fig.  4. 


Consider  a stream  of  gas  incident  on  the  plates  coming  from  the 


direction  defined  by  azimuthal  and  polar  angles  cpv  and  0 v.  The  unit 


vector  in  the  direction  of  cpv  and  0V  is  given  by 


sin  0 cos  m x + sin  G„  cos  cp  y + cos  0V  £ 

V V 1 V v ~ v 


Thus  the  velocity  of  the  gas  molecules  is  - v oo  I v.  The  magnitudes  of 


the  drag  and  lift  forces  on  plate  1 are  then  given  ty  Equations  (4)  and  (5) 


with  cos  y = v . n.  . If  cos  y < 0,  the  gas  stream  is  incident  on  the  back 


side  of  the  plate,  in  which  case  both  drag  and  lift  forces  will  be  considered 


to  be  zero.  (The  plate  is  shadowed  by  other  surfaces  of  the  satellite.  ) 


^ ^ ^ ^ ^ 

Unit  vectors  in  the  drag  and  lift  directions  are  Fp  = -v  and  Fj^  = (njX  v)  x v. 


Denoting  either  drag  or  lift  vectors  by  F (^),  a spin  averaged  force  acting 


on  the  pair  of  plates  is  given  by 


:F>  = — 


* / 


?o  + IL 

r 2 


F (?)  d?  , 


So  - E 


A A 

where  ^Q=  cp^is  the  spin  angle  for  which  n j is  closest  to  v. 
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Equation  (9)  was  used  to  calculate  spin-averaged  dr&g  and  lift  forces 
for  the  S'3-l  satellite.  The  asymmetric  shape  due  to  the  inclined  "roof" 
causes  a non-zero  average  lift  force.  The  ratio  of  lift  to  drag  forces 
for  various  values  of  a and  P is  shown  in  Fig.  5.  It  will  be  shown  later 
that  only  values  of  a near  1 appear  tobe  reasonable.  The  dependence  on  P 
of  the  lift  to  drag  force  ratio  for  large  values  of  a is  seen  to  be  weak. 

In  addition,  other  investigations  of  satellite  drag  phenomena  have  given 
results  compatible  with  diffuse  rather  than  specular  reflection  (Imbro, 
et  al.,1975;  Beletskii,  19701.  Hereafter  we  will  therefore  consider  only,  the 
case  P = 1 (diffuse  reflection). 

The  nominal  orientation  of  the  spin  axis  perpendicular  to  the  orbit 
plane  was  maintained  to  within  about  5 degrees.  To  test  for  sensitivity 
to  departures  from  nominal  orientation,  the  ratio  of  spin  averaged  lift  to 
drag  force  was  calculated  fora  range  of  pitch  angles  on  either  side  of 
the  nominal.  The  results,  shown  in  Fig.  6,  indicate  a strong  dependence 
on  actual  spin  axis  orientation,  a pitch  angle  of  only  4 deg  being  sufficient  to 
reduce  the  lift  force  to  zero.  For  the  calculations  described  below  we 
have  therefore  used  the  flight -determined  spin  axis  location  rather  than 
assuming  nominal  orientation.  Only  the  availability  of  frequently  determined, 
high  quality  attitude  and  ephemeris  data  has  allowed  reliable  interpretation 
of  variations  of  orbital  elements,  for  this  satellite  with  its  large  lift  to  drag 
force  ratio. 
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RATIO  LIFT/DRAG,  PERCENT 


EFFECTS  OF  WINDS  AND  LIFT  ON  ORBITAL  INCLINATION 


To  determine  atmospheric  winds  from  changes  in  orbital  inclination, 
it  was  first  necessary  to  correct  the  data  for  other  known  effects. 
Specifically,  lunisolar  and  geopotential  zonal  harmonic  perturbations  were 
removed;  solar  radiation  pressure  was  found  to  be  negligible,  and  no 
corrections  were  made  for  it. 

In  principle,  aerodynamic  lift  effects  could  have  been  removed  from 
the  inclination  data  in  the  same  manner  that  the  other  perturbations 
were  handled.  It  was  found,  however,  that  the  magnitude  of  the  lift 
effect  was  very  sensitive  to  the  value  adopted  for  the  accommodation  co- 
efficient, which  is  not  well  known  or  understood.  It  was  decided,  there- 
fore, that  the  wind  analysis  needed  to  be  conducted  for  a range  of 
reasonable  values  of  the  accomodation  coefficient.  The  derivation  of  the 
lift  force  was  presented  in  the  previous  section.  Here,  the  lift  force  is 
assumed  to  be  known,  and  its  effect  on  the  inclination  will  be  discussed. 

A particular  problem  that  existed  in  the  present  study  was  that  the 
orbital  eccentricity  in  the  early  portion  of  the  data  span  exceeded  a 
value  of  0.20.  This  precluded  the  straightforward  application  of  King- 
Hele's  well-known  analytic  expressions  (for  small  eccentricities)  that 
relate  changes  in  orbital  inclination  to  atmospheric  winds  (e.g.  King- 
Hele,  1966).  Similarly,  in  deriving  the  inclination  change  due  to 
aerodynamic  lift,  the  "small  eccentricity"  approximation  was  avoided. 


The  problem  of  "large  eccentricities"  has  recently  been  treated 
analytically  by  King-Hele  and  Walker  (1976),  but  as  their  work  was  not 
available  to  us  at  the  time  of  our  analysis,  we  shall  proceed  to  describe 
our  approach,  which  was  numerical  rather  than  analytic. 

The  theory  relating  orbital  inclination  changes  to  atmospheric  winds 
and  other  forces  normal  to  the  orbit  has  been  presented  in  detail  by 
King-Hele  (e.  g.  1964,  1966);  therefore  no  attempt  will  be  made  here  to 
reproduce  his  thorough  work.  It  shall  be  assumed  that  the  reader  is 
familiar  with  King-Hele's  work,  and  only  major  points  of  relevance  to 
the  present  study  will  be  pursued. 

Changes  in  the  orbital  inclination  arise  only  from  forces  normal  to 
the  orbit  plane.  In  an  orbital  period  T,  the  average  change  Ai  resulting 
from  a force  per  unit  mass  f that  is  normal  to  the  orbit  may  be 
written  (King-Hele,  1964) 


A 


l 


T 

r f cos 


a (1 


u dt 


(10) 


where  r is  the  geocentric  distance  of  the  satellite,  p is  the  gravita- 
tional constant  of  the  earth,  a is  the  semi-major  axis,  and  e is  the 
eccentricity.  The  angle  u is  the  argument  of  latitude,  equal  to  the  sum 
of  the  argument  of  perigee  m and  the  true  anomaly  0.  The  forces 
normal  to  the  orbit  that  are  of  concern  are  those  due  to  aerodynamic 
drag  (i.e.  atmospheric  winds)  and  aerodynamic  lift. 
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Drag.  Using  the  eccentric  anomaly  E instead  of  time  tas  the 
independent  variable,  King-Hele  (1966,  equation  35)  has  shown  that  the 
change  in  inclination  per  orbital  revolution  arising  from  the  drag  force 
normal  to  the  orbit  is 


rZn  2 r - 1/2 

I r exp(z  cos  E - z)  r(l  + e cos  E) 
cT  vp(l  - e2) 

a sin  i cos2  (uj  + 0)  - $ cos  i cos  + 0) 

£l  - sin2  i sin2(g)  + 6)J^ 


where  Dp  is  the  density  at  perigee,  and  z = ae/H,  H being  the  density  scale 
height  at  perigee.  (Because  ot  the  high  eccentricity  of  the  orbit,  constant 
scale  height  and  spherical  atmosphere  assumptions  are  utilized.  ) The 
quantities  w and  $ are,  respectively,  the  zonal  and  meridional  angular 
velocities  of  the  atmosphere,  measured  positive  eastward  and  north- 
ward. The  function  y is  nearly  constant  and  may  be  taken  to  be 

yv”  3 1 - rp  w cos  i / vp , 

where  r^  js  the  perigee  distance  and  v^  the  velocity  at  perigee.  The 
quantity  6 in  Equation  (11)  is  a function  of  the  physical  characteristics 
of  the  satellite: 


6 = y S C /m 

D D s, 


1 


where  C^is  the  drag  coefficient  appropriate  for  a reference  area 

of  the  satellite,  and  m is  the  mass  of  the  satellite. 

s 

The  approach  used  by  King-Hele  (1966)  to  evaluate  the  integral  in 
Equation  (11)  involved  expansions  of  the  functions  in  powers  of  e,  thereby 
limiting  the  resulting  expression  for  <Ji  to  small  values  of  e.  While 
this  approximation  is  adequate  for  the  great  majority  of  satellites,  the 
error  would  be  too  large  for  this  particular  satellite,  and  thus  numerical 
integration  was  employed  in  this  study. 

As  pointed  out  by  King-Hele,  it  is  useful  to  eliminate  the  explicit  de- 
pendence of  Ai  (Equation  11)  on  the  perigee  density  by  use  of  the  change 
in  orbital  period  per  revolution,  AT.  It  may  be  shown  (King-Hele,  1964) 


AT  = - 6tt2  7(a/p)  6 a2  Pt 


e 3 B(z) 


where  B(z)  = IQ  + 2elf  + 3 (IQ  + I^)  + e ^ (3  + I. 


The  I ' s in  the  expression  above  are  the  modified  Bessel  functions  of 
the  first  kind  of  argument  z a e/H.  It  should  be  noted  that  Equation  (12) 
was  also  derived  by  expansion  in  powers  of  e,  but  that  sufficient  higher 
order  terms  were  retained  to  permit  its  use  without  serious  error  for 
eccentricities  e ~ 0.  2. 
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Substitution  of  Equation  (12)  into  (11)  leads  to 


6ttB(z) 


Ztt 

/«>■ 


r (1  + e cos  E) 
a (1  - e 2) 


1/2 


exp(z  cos  E)  x 


A sin  i cos2  (uy  + 0)  - e cos  i cos(w  + 8)  dE, 

j^l-sin2  i sin^  (m+G)  * ^ 


(13) 


where  ATyis  the  change  in  orbital  period  per  revolution  expressed  as 
a fraction  of  a day,  and  A and  e are  related  to  the  zonal  and  meridional 
wind  components  by 

A - w / wg 
e - 0/we  , 

where  we=  2 it  radians  per  day.  It  is  to  be  noted  that  the  ratio  r/a 
and  the  true  anomaly  0 are  functions  of  E: 

r/a  = 1-  e cos  E , 

sin  0 = (a/r)  \A-e2  sin  E , 

cos  0 = (a/r)  (cos  E - e) 

Lift.  Because  of  attitude  control,  the  spin-averaged  lift  force  on  this 
satellite  was  at  all  times  within  a few  degrees  of  the  normal  to  the  orbit. 
The  magnitude  of  the  lift  force,  like  that  of  the  drag  force,  is  proportional 
to  the  square  of  the  satellite  velocity  relative  to  the  ambient  air. 
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However,  little  error  results  if  the  magnitude  is  taken  to  be  proportional 
to  the  square  of  the  inertial  velocity.  The  force  per  unit  mass  normal 
to  the  orbit  produced  by  lift  may  then  be  written 


f 


L = 


1 

2 


pv2(^  Cl  /ms) 


(14) 


where  p is  the  atmospheric  density,  v is  the  inertial  velocity  of  the 
satellite,  and  is  the  coefficient  of  lift  appropriate  for  an  effective 
satellite  area  Sj^  . 

By  use  of  Equations  (10)  and  (14)  and  the  relationship 
2 

v = n (1  + e cos  E)  / r. 


as  well  as  other  relationships  between  the  orbital  parameters,  it  may 
be  shown  that  the  change  in  inclination  per  revolution  due  to  lift  is 
given  by 

Ai  = - tt  p a(ST  C /m  )e’ZA  (z)(l-e2)  2 cos  ou,  (15) 

L P L L s 


whe  re 


A(z ) = (1-e'  ) I 


, - e T + e I, 

1 2 0 2 2 


and,  as  before,  the  I are  the  modified  Bessel  functions  of  the  first 

n 

kind  of  argument  z. 
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Using  Equations  (10),  (12)  and  (14)  and  letting  be  the  orbital 
period  in  units  of  a day  one  obtains 

Ai^  = 2VTd  cos  uj  Sl  Cl  A(z)  (16) 

3td  /i 

In  Equation  (16),  deviation  of  the  satellite  spin  vector  from  the 
nominal  position,  which  would  affect  the  lift  force  as  described  in  the 
previous  section,  is  taken  into  account  by  the  term  ^as 

been  assumed  constant  over  an  orbital  period,  but  varies  on  a longer  time 
scale. 

Finally,  in  more  concise  form,  the  sum  of  drag  (i.  e.  atmospheric 
wind)  and  lift  effects  on  the  orbital  inclination  is  expressed  by 

Ai  = Y (of)  + AGz  + e Gm>  (17) 

where  Y (a)  is  the  aerodynamic  lift  contribution,  given  by  AiL  /ATp 

((Equation  16),  whose  dependence  on  a,  the  thermal  "accommodation 

coefficient,"  arises  from  the  ratio  . The  ratio  of  lift  to 

drag  forces,  discussed  in  the  previous  section,  yields  (S^C-j^/S^Cp).  The 
second  and  third  terms  on  the  right  hand  side  of  Equation  (17)  represent 
the  7. onal  and  meridional  wind  contributions,  respectively;  expressions  for 

C.  and  G are  readily  derived  from  Equation  (13). 

7.  m 
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RESULTS 


Equation  (17)  relates  the  rate  of  change  of  inclination  to  the  orbital 
parameters,  the  wind  speeds  (Aand  e),  the  thermal  accommodation 
coefficient  a and  the  satellite  geometry  and  orientation.  Since  all  of 


these  quantities  except  A,  e and  a are  known,  values  of  A and  e may  be 
determined  which  "best"  fit  the  data  corrected  for  luni-solar  and  geo- 
potential perturbations  for  a given  value  of  a.  This  best  fit  was  obtained 
by  integrating  Equation  (17)  and  evaluating  Aand  e and  the  initial  value 
of  i in  such  a way  that  the  data  were  fit  in  a least- square s sense.  A total 
of  213  data  points  were  analyzed  in  this  manner,  and  the  results  of  this 
procedure  are  described  in  this  section.  As  shown  by  King-Hele,  the 
inferred  atmospheric  wind  speeds  apply  to  an  altutude  3/4  of  a scale 
height  above  the  perigee  altitude  of  ~l60km. 

In  this  analysis,  the  data  set  was  treated  in  two  ways.  First, 
least  - squares  fits  v/ere  carried  out  for  the  entire  data  set  for  six  values 
of  C : 1.  0.  0.  99,  0.97,  0.  95,  0.85  and  0.70.  (Note  that  the  effect  of 
aerodynamic  lift  vanishes  for  O'  1.0.)  In  all  cases,  it  was  found  that 
the  value  of  e was  poorly  determined,  and  variations  in  the  value  of  e 
had  little  effect  on  the  value  of  A deduced.  This  is  primarily  due  to  the 
fact  that  i is  near  90  deg  so  that  the  effect  of  e on  i is  small  (cf.  Equation 
13).  Therefore,  the  fits  were  made  with  e = 0,  and  no  further  discussion 
of  is  given  here.  Figure  7 shows  the  least- squares  fit  to  the  data  for 

>j  1.0,  for  which  the  best  fit  value  of  A is  1.  31  (solid  curve  ) and  for  Of 

0,99,  for  which  the  best  fit  value  of  A is  1.22  (dashed  curve).  Both  cal- 
culated curves  reproduce  the  general  evolution  of  the  inclination  with  time 
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the  greatest  rate  of  change  of  inclination  occurring  during  periods  in  which 
the  satellite  perigee  was  near  the  earth's  equator  (uu  =0,  180  deg)  and 
near  the  end  of  the  satellite's  lifetime,  when  the  orbital  eccentricity  de- 
creased rapidly.  However,  quantitative  comparison  indicates  that  the  cal- 
culated best  fits  depart  substantially  from  the  data  (i.  e.  , differences  on 
the  order  of  0.01  - 0.03  deg). 

Table  1 gives  the  calculated  values  of  A for  the  six  values  of  a noted 
above,  along  with  the  rms  error  of  the  fit.  Several  features  are  evident. 
The  inferred  value  of  A is  fairly  sensitive  to  changes  in  a,  especially  for 
values  of  a very  near  1.  In  fact,  the  value  of  A inferred  from  this  analysis 
indicates  an  average  eastward  wind  relative  to  the  rotating  earth  if 
a>0.  85  and  a westward  wind  otherwise.  However,  inspection  of  the  value 
of  the  rms  error  as  a function  of  <y  discloses  that  the  rms  error  minimizes 
in  the  range  0.95  ' a < 0.  99,  so  that  values  of  a in  this  range  are  more 
acceptable  than  others.  Orbital  considerations,  coupled  with  theortical 
treatments  of  thermospheric  winds  (Blum  and  Harris,  1975;  Straus  et  al.  , 
1975a,  b),  also  tend  to  favor  the  value  of  A as sociated  with  the  larger  values 
of  a.  As  noted  earlier,  the  local  time  of  satellite  perigee  lies  nearly  in 
the  8 hr  local  time  - 20  hr  local  time  plane  (see  Table  2).  During  the  life- 
time of  the  satellite,  perigee  occurred  on  the  evening  side  of  the  earth 
approximately  67%  of  the  time,  so  that  the  overall  value  of  i is  weighted 
twice  as  strongly  by  its  value  at  ~20  hr  local  time,  where  the  theoretical 
studies  cited  above  indicate  a strong  eastward  wind,  than  by  that 


Table  1. 


Calculated  values  of  A for  the  entire  data 
set  for  several  assumed  values  of  <y 


RMS 

error  (deg) 

1.  24 

X 

10-2 

7.  87 

X 

10-3 

6.  59 

X 

>— < 
o 

1 

UJ 

7.  00 

X 

10-3 

1.  19 

X 

10'2 

1.74 

X 

10‘2 

1.  12 
1.  00 
0.  91 


Table  2.  Characteristics  of  the  six 
data  segments 


SEGMENT 

HEMISPHERE 

LST 

NO. 

POINTS 

1 

Northern 

21-20 

60 

2 

Southern 

20-21 

38 

3 

Southern 

8 

37 

4 

Northern 

7 

32 

5 

Northern 

19 

31 

6 

Southern 

19 

15 

at  ~8  hr  local  time,  where  westward  winds  are  expected.  Thus,  the 


calculated  value  of  A~l.10-1.20  is  understandable,  and  does  not  neces- 
sarily indicate  an  average  eastward  motion  of  the  atmosphere. 

The  excellent  quality  of  the  orbital  data  allows  more  detailed 
analysis  of  the  changes  in  inclination  than  that  described  above.  In  order 
to  investigate  diurnal  variations  of  A,  and  to  further  elucidate  the  effects 
of  different  assumed  values  of  a on  the  inferred  values  of  A,  the  data  set 
has  been  treated  in  a second  manner,  as  follows.  As  shown  in  Table  2 
and  in  Fig.  8,  the  data  set  was  divided  into  6 segments  in  each  of  which 
the  perigee  latitude  was  restricted  to  lie  in  one  hemisphere  and  also  to 
occur  on  one  side  of  the  earth.  For  example,  Segment  1 contains  data 
points  for  which  perigee  occurred  in  the  northern  hemisphere  on  the 
evening  ( ~22  hr  local  time)  side  of  the  earth,  whereas  Segment  2 

covers  the  southern  hemisphere  on  the  evening  side.  This  division  thus 
allows  us  to  distinguish  variations  of  A with  local  time  as  well  as  season. 

The  least- squares  fit  procedure  was  carried  out  for  each  of  the  six 

segments  of  data  for  four  values  of  a : 1.  0,  0.  99,  0.  95  and  0.85.  An 

example  of  the  accuracy  of  the  fit  is  given  in  Fig.  8,  which  displays 

every  second  data  point  as  an  "x",  the  calculated  fit  for  a - 0.99  as  a 

piecewise  continuous  solid  line,  and  the  residual  for  each  data  point.  The 

-3 

rms  error  for  the  entire  data  set  in  this  case  is  1.  33  x 10  deg,  a factor 
of  almost  6 better  than  that  obtained  when  the  entire  data  set  was  treated 
at  once.  This  overall  agreement  is  considerably  better  than  that  which 
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has  been  attained  in  this  type  of  analysis  in  the  past. 

The  calculated  values  of  A for  each  of  the  six  segments  are  given  in 
Table  3.  In  each  segment,  a substantial  variation  with  a of  the  inferred 
value  of  A is  evident,  and,  for  a fixed  value  of  a,  the  computed  value  of 
Avaries  from  segment  to  segment.  There  is  no  a_  priori  way  to  pick  the 
correct  value  of  a.  Examination  of  the  residuals  of  the  fit  shows  no  strong 
tendency  towards  a particular  value  of  a,  although  the  fits  with  o?>  0.  95 
are  consistently  better  than  those  for  a - 0.85;  however,  in  no  case  were 
the  residuals  for  a 1.0  less  than  those  for  a - 0.99.  On  the  other  hand, 

if  a single  value  of  a applies  for  the  entire  lifetime,  the  values  of  A 
inferred  for  cr>  0.  99  are  preferable  on  the  basis  of  comparision  with  the 
theoretical  studies  cited  above.  These  studies  predict  an  eastward  (A>1) 
motion  durirg  the  evening  (Segments  1,2,5  and  6)  and  a westward  (A<  1 ) 
motion  during  the  morning  (Segments  3 and  4);  only  the  results  for  values 
of  q very  near  1.0  show  this  behavior. 
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DISCUSSION  AND  SUMMARY 


The  present  paper  has  dealt  with  the  effects  of  upper  atmospheric  winds 
and  aerodynamic  lift  on  the  changes  in  inclination  of  the  orbit  of  the  S3-1 
satellite.  Aerodynamic  lift  effects  have  not  been  considered  in  previous 
studies  to  determine  atmospheric  winds  from  changes  in  inclination,  and 
only  the  high  quality  of  orbit  and  attitude  data  for  this  satellite  has  per- 
mitted it  in  this  case. 

The  argument  is  occasionally  advanced  that  aerodynamic  lift  effects 
may  be  ignored  for  satellites  which  do  not  have  an  active  orientation 
mechanism  because  such  satellites  undergo  random  tumbling  motion.  It 
should  be  noted,  however,  that  in  the  absence  of  disturbances  such  as 
aerodynamic  torque,  a spinning  satellite  will  maintain  its  angular  momentum, 
and  thus  its  spin  axis  will  be  inertially  fixed.  For  asymmetric  satellites 
(such  as  S3-1)  the  spin  averaged  lift  force  may  be  non-zero.  Great  care 
should  therefore  be  exercised  before  aerodynamic  lift  effects  are  ignored 
for  asymmetric  satellite  shapes,  even  for  those  which  are  tumbling 
"randomly.  " 

As  discussed  in  the  previous  section,  the  lift  force,  and  hence  the 
inferred  values  of  the  zonal  wind  speed  and  direction,  are  quite  sensitive 
to  the  assumed  value  of  the  thermal  accommodation  coefficient.  Analysis 
of  the  zonal  winds  derived  using  different  assumptions  of  the  value  of  a 
indicates  that  best  agreement  with  current  theoretical  thermospheric 
models  is  attained  when  nearly  complete  thermal  accommodation  is 
assumed  (i.  e,  a « 1.  0).  There  have  been  a few  studies  dealing  with  the 
value  of  a appropriate  for  an  orbiting  satellite.  Imbro  et  al.  (1975), 
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using  spin  rate  decay  data  from  Ariel  2(perigee  altitude  ~300  km),  estimated 
an  accommodation  coefficient  near  0.  9 and  discussed  the  possible  increase 
of  a with  decreasing  altitude  (associated  with  the  increase  in  atmospheric 
concentration  of  atomic  oxygenj.  In  fact,  Beletskii  (1970)  has  analyzed 
spin  rate  data  from  the  paddle  wheel  satellite  Proton  2 (perigee  altitude  ~ 

190  km)  and  concluded  that  a > 0.  99.  Beletskii' s results  and  the  discussion 
of  Imbro  et  al.  thus  support  the  assumption  of  large  values  of  o,  which 
seem  most  reasonable  on  the  basis  of  the  effects  of  aerodynamic  lift  on 
inferred  atmospheric  wind.  These  values  of  a were  derived  assuming 
models  which  did  not  include  possible  surface  chemistry  effects.  Atomic 
oxygen,  a major  constituent  in  the  thermosphere,  is  highly  reactive,  and 
it  is  not  unlikely  that  a satellite  surface  may  become  covered  with  a chemi- 
sorbed layer  of  atomic  oxygen.  An  incoming  atomic  oxygen  atom  might 
then  collide  with  the  surface,  recombine  with  a surface  adatom,  and  be 
re-emitted  as  a molecule.  Present  knowledge  of  surface  chemistry  does 
not  allow  treatment  of  such  effects  with  confidence. 

Assuming  that  values  of  a » 1 are  indeed  correct,  the  present  study 
provides  some  information  regarding  the  diurnal  variations  of  the  zonal 
wind  speed  and  direction  at  ~175  km  altitude.  For  local  times  in  the 
range  19-21  hr  the  calculated  wind  is  eastward  and  has  a magnitude  of 
~ 400  m/sec.  For  local  times  of  7-8  hr,  the  calculated  wind  is  west- 
ward and  on  the  order  of  200  m/sec.  These  directions  are  in  agreement 
with  the  theoretical  studies  cited  earlier,  and  the  wind  speeds  are  in 
general  agreement  with,  or  somewhat  larger  than,  those  predicted  theo- 
retically. The  wind  speeds  also  are  similar  to  the  results  of  other  obser- 
vations at  comparable  altitudes  (Ching,  1971;  King-Hele,  1972b;  Forbes, 

1975). 
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No  seasonal  variations  in  wind  speed  are  evident  in  the  present  results  with 


O'  ~ 1.  The  possible  dependence  of  wind  speed  on  the  level  of  magnetic 
activity,  noted  by  Forbes  (1975),  could  not  be  investigated,  since  the 
average  value  of  Ap  for  each  of  the  6 data  segments  was  about  the  same. 

Since  the  value  of  o is  not  well  known,  and  since  the  inferred  wind 
speeds  and  directions  are  so  sensitive  to  small  deviations  of  a from  unity, 
the  above-mentioned  wind  speeds  can  be  regarded  as  being  only  approximate. 
Although  it  is  gratifying  that  "acceptable"  wind  properties  have  been  found 
to  be  consistent  with  "acceptable"  accommodation  coefficients,  it  is 
nevertheless  apparent  that  accurate  atmospheric  wind  speeds  cannot  be 
obtained  from  this  particular  satellite  because  of  the  uncertainty  in  the 
accommodation  coefficient.  The  effects  of  lift  and  its  dependence  on  the 
accommodation  coefficient  will,  of  course,  depend  on  the  individual 
characteristics  of  a satellite's  shape  and  orientation.  The  almost  wing- 
like shape  of  S3-1  undoubtedly  is  capable  of  generating  lift  forces  larger 
than  would  more  symmetrically  shaped  satellites.  If  accurate  zonal  wind 
velocities  are  to  be  derived  from  measurements  of  satellite  orbital  incli- 
nation changes,  care  must  be  taken  to  ensure  that  the  shape  and  orientation 
of  the  satellite  are  such  that  the  lift  forces  generated  are  sufficiently  small 
that  the  resulting  uncertainty  due  to  inadequate  knowledge  of  the  accommo- 
dation coefficient  does  not  cause  significant  error. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  acientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment ; application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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